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Introduction 

The fractional quantum Hall effect (FQHE) occurs in a two-dimensional elec- 
tron gas (2DEG) in a strong magnetic field oriented perpendicular to the plane of 
the electronsfl]. The effect was discovered as a transport anomaly. In a transport 
measurement it is noted that at certain strengths B*(n), which depend on the density 
n of the 2DEG, current can flow without any dissipation. That is, there is no voltage 
drop along the flow of the current. At the same time, the Hall voltage perpendicular 
to both the direction of the current and of the magnetic field is observed to attain a 
quantized value for a small, but finite, range of magnetic field or density, depending on 
which quantity is varied in the experiment. The effect is understood to be the result 
of an excitation gap in the spectrum of an infinite 2DEG at these magnetic fields. A 
convenient measure of the density of a 2DEG in a strong magnetic field is given by 
the filling factor v = 2ir£^n, with £b = Jhc/(eB) the magnetic length. The filling 



factor gives the ratio of the number of particles to the number of available states in 
a magnetc sub-band (Landau level), or, equivalently, the number of particles per flux 
quantum $o — hc/e. The quantum Hall effect was first discovered[2] at integer filling 
factors. In this integer quantum Hall effect, the energy gap is nothing but the kinetic 
energy gap hu c = HeB/(m*c). Later, the fractional quantum Hall effect (FQHE) was 
discovered[3] at certain rational filling factors of the form v = p/q, with p and q relative 
primes, and q odd. In the FQHE, the excitation gap is a consequence of the strong 
electron-electron correlations. Therefore, any computational approach to the quantum 
Hall effect must accurately treat the electron correlations in order to capture the FQHE 
at all. We emphasize that this fact makes any density functional approach to the FQHE 
qualitatively different from 'standard' applications to atomic, molecular or condensed 
matter systems. In such standard applications, the electron correlations usually give 
a quantitatively important correction of the order of perhaps 10 - 20% but do not 
usually qualitatively change the results in a fundamental way. Furthermore, standard 
applications of density functional theory do not require ensemble density functional 
theory (in such applications, the systems under considerations are typically pure-state 
f-representable). Our work represents (to the best of our knowledge) the first practical 
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applications of ensemble density functional theory to a system which is not pure-state 
f-representable. 

Our understanding of the origin of the FQHE started with Laughlin's seminal paper 
of 1983 [4], which dealt with the simplest fractions v = 1/m, with m an odd integer. 
At these values of u, there are on the average m magnetic flux quanta <£>o = hc/e per 
electron. Laughlin constructed a variational wavefunction for spin-polarized systems 
in strong magnetic fields, strong enough that the splitting hu c between the magnetic 
subbands, or Landau levels, can be taken to be infinite. The wavefunction can then be 
constructed from single-particle states entirely within the lowest Landau level. Laughlin 
wrote the variational wavefunction as 



where Zj = Xj + iyj is the coordinate of the jth electron in complex notation. This 
wavefunction is an eigenstate of angular momentum. Laughlin demonstrated that the 
system having the wavefunction Eq. ([]]) is an incompressible liquid with v = 1/m, m 
odd, with an energy gap to excitations, and that the elementary excitations are frac- 
tionally charged quasi- holes or quasi-particles of charge e* = ±e/m. The origin of the 
energy gap can be understood in the so-called pseudo-potential representation of the 
electron-electron interactions [5]. Here, the electron-electron interaction V(r$ — r^) be- 
tween electrons i and j is decomposed into strengths Vi in relative angular momentum 
channels i = 0, 1, 2, . . ., of the two electrons. For any realistic interaction V(ri — rj), it 

turns out that Vq > Vx > Vi > Consider the case v = 1/3. In this case, the lowest 

angular momentum pseudo-potential that enters into the m = 3 Laughlin description 
is V\, the interaction energy of two electrons of unit relative angular momentum (in 
units of fi). (Only odd relative angular momenta are permissible for spin-polarized 
electron wavef unctions, since they have to be anti-symmetric under interchange of elec- 
tron coordinates.) The Laughlin wavefunction is a very cleverly constructed highly 
correlated state which completely excludes unit relative angular momentum between 
any two electrons, and is furthermore the only state which satisfies this property at 
v — 1/3. Therefore, any excited state at this filling factor must contain some electrons 
with unit relative angular momentum. The energy gap is due to the cost of this, and 
hence is of order V\. Figure [I] depicts the exchange-correlation energy per particle for 
infinite, homogeneous FQHE systems vs. filling factor. The cusps at filling factors 
v = 1/3,2/5,3/5,2/3,4/3, and v = 7/5 have been included to scale. Note that these 
cusps are barely visible on this scale, yet they are responsible for all the physics of the 
FQHE! 

Finite systems, and spins 

We have outlined above how the electron-electron interactions in an infinite, ho- 
mogeneous system produce the excitation gap. It is important to note that these gaps 
are only for excitations in the bulk of the system. When a system is bounded there 
must be gap less excitations located at the boundaries of the system [6, 7]. Since all 
experimental systems are finite and inhomogeneous, the low-energy properties probed 
by experiments are necessarily determined by the gapless edge excitations. Advances in 
semiconductor nanofabrication technologies have lead to the possibility of manufactur- 
ing systems which are extremely inhomogeneous, and in practice dominated by edges. 
As an example, recent experiments have even been performed on tiny quantum dots, 
containing about 30 electrons[8, 9]. Also, edge structures in inhomogeneous systems 





2 



0.0 




-0.8 1 1 1 1 1 1 1 

0.0 0.5 1.0 1.5 

Filling factor v(r) 

Figure 1. The ground state energy e xc per particle of an infinite, homogeneous, spin-polarized 
FQHE system is depicted as a function of filling factor. The cusps at v = 1/3, 2/5, 3/5, 2/3, 4/3, and 
v = 7/5 are included. In this plot, the smooth part of the exchange-correlation energy was obtained 
from the Levesque-Weiss-MacDonald interpolation formula Eq. ([?]) 

have been studied directly using capacitance spectroscopy [10], time-resolved measure- 
ments of edge magnetoplasmons (the gapless modes along the boundary) [11], and by 
surface acoustic wave techniques capable of resolving very small inhomogeneities in the 
electron density [12]. 

So far, we have limited the discussion to spin-polarized system. At first look it 
seems reasonable that in the strong magnetic fields used in experiments on QHE sys- 
tems the Zeeman splitting g\i B B, where g is the effective Lande factor and \xb the Bohr 
magneton, is large enough that the high-energy spin direction is energetically inacces- 
sible. However, two factors conspire to make the Zeeman splitting very low in GaAs 
systems. First of all, spin-orbit coupling in the GaAs conduction band effectively lowers 
the Lande factor to g ~ 0.44. Second, the low effective mass, m* « 0.067m e , further 
reduces the ratio of Zeeman energy to cyclotron energy to about 0.02, compared to its 
value of unity for free electrons. For magnetic fields of about 1 - 10 T, the Coulomb 
energy scale of the electron-electron interactions, e 2 /(e ^), where e ~ 12.4 is the 
static dielectric constant, is of the same order as the cyclotron energy. As a first ap- 
proximation, one should then set the Zeeman energy to zero, rather than infinite, since 
it is two orders of magnitude smaller than the other energy scales. As a consequence, 
the spin-degree of freedom is governed by the electron-electron interactions, rather 
than by the Zeeman energy. This dramatically changes the nature of the low-energy 
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bulk single-particle excitations near filling factors of the form v = 1/m, with m odd, 
from single-particle spin-flips to charge-spin textures. In these objects, loosely called 
skyrmions, the spin density varies smoothly over a distance of several magnetic lengths, 
so that the system can locally take advantage of the exchange energy. Skyrmions in the 
quantum Hall effect were first predicted by Sondhi et al. [13]. Recent experiments [14] 
measuring the spin-polarization near v — 1 have provided evidence for this. According 
to calculations by Karlhere et al. [15], the skyrmion excitations in an infinite homoge- 
neous system may lead to charge-spin textured appearing at the edges of QHE systems, 
which would dramatically change our understanding of QHE edges. 

In order to accurately understand the experiments and inhomogeneous FQHE sys- 
tems in general, we must have a way of accurately calculating their properties, including 
the spin degrees of freedom. Certain aspects of inhomogeneous FQHE systems have 
been studied by different techniques. For example, field theories can be constructed to 
study the low-energy limit of the gapless edge excitations for spin-polarized systems [16], 
or the spin texture on the edge of an infinite Hall bar [15]. However, effective field the- 
ories typically do not give accurate quantitative results, other than for e.g. critical 
exponents near fixed points. So-called composite fermion methods have been used for 
non-interacting composite fermions[17] and in a Hartree approximation to study finite 
spin-polarized FQHE systems[18, 19, 20]. In this approach, the Chern-Simons term, 
arising from the singular gauge transformation, is replaced by its smooth spatial aver- 
age, and the composite fermion mass, which is not well known, has to be put in by hand. 
The electron-electron interaction can then treated in a self-consistent Hartree approxi- 
mation. The hope is then that the most important aspects of the electron-electron corre- 
lations are included in this approximation. Near v — 1, at which the Slater determinant 
is the exact ground state, it makes sense to use the Hartree-Fock approximation, 
and the stability of a spin-polarized quantum dot at v = 1 as a function of confining 
potential has been studied in this approximation [21, 22]. The edge structure of spin- 
polarized systems has also been studied using semiclassical methods[23, 24], in which 
the electron-electron interaction is included at the Hartree level and it is furthermore 
assumed that all potentials vary on a length scale much larger than £ B . Beenakker[23], 
and Chklovskii, Shklovskii and Glazman[24] demonstrated that the edge of an integer 
quantum Hall system, in which the correlation energies between the electrons can be 
ignored, consists of a sequence of compressible and incompressible strips. The widths 
of the incompressible strips is determined by the length over which the effective con- 
fining potential (external plus Hartree) varies an amount equal to the energy gap hu c . 
The origin of the compressible and incompressible strips are the energy gaps, which are 
the kinetic energy gaps hou c in the case of the integer quantum Hall effect. But it is 
easy to generalize the argument to include the energy gaps causing the FQHE[24]. The 
conclusion is then that there should be compressible and incompressible strips, with 
the density of each incompressible strip fixed at the value of an FQHE fraction. The 
width of each incompressible strip is then fixed by the length over which the effective 
confining potential varies an amount equal to the energy gap of the FQHE fraction 
corresponding to the density of that strip. 

Quantum dots including the spin degree of freedom have also been studied by 
direct numerical diagonalizations[25]. These calculations demonstrated the importance 
of the spin degree of freedom. At the present, numerical diagonalizations are limited 
to systems with of the order of 10 electrons. 

It is highly desirable to have a computational approach which accurately includes 
electron-electron correlations and spin degree of freedoms, and which can handle inho- 
mogeneous systems with on the order of 10 2 -10 3 electrons. One such approach which 
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is in principle valid for any interacting electron system is density functional theory 
(DFT)[26, 27, 28]. There have been some attempts to apply density functional theory 
to the FQHE. Ferconi and Vignale[29] applied current density functional theory[30] to 
small, parabolically confined quantum Hall systems and showed that the current density 
functional theory gave good results for the ground state energy and spin polarization 
near v — 1. However, the energy gaps due to correlation effects were not included in 
that calculation. Ferconi, Geller and Vignale[31] also recently studied spin-polarized 
FQHE systems within the spirit of the DFT using an extended Thomas- Fermi approxi- 
mation at low, but non-zero, temperatures. In this, the kinetic energy was treated as a 
local functional, as in the standard Thomas-Fermi approximation, while the exchange- 
correlation energy was included in a local density approximation (LDA). This extended 
Thomas-Fermi approximation is valid in the limit of very slowly varying confining po- 
tential. Ferconi, Geller and Vignale focused on the incompressible and compressible 
strips at an edge of an FQHE system, and obtained results in agreement with the 
predictions by Chklovskii, Shklovskii, and Glazman[24]. 

We have developed for the fractional quantum Hall effect an ensemble DFT scheme 
within the local density approximation, and have applied it to circularly symmetric 
quantum dots [32, 33]. In our approach, the kinetic energy is treated exactly, and the 
density represented by Kohn-Sham orbitals. The results are in good agreement with 
results obtained by semiclassical[23, 24, 31], Hartree-Fock[21, 22] (for cases where the 
correlations do not play a major role), and exact diagonalization methods [25]. Our 
calculations for spin-polarized systems show that the exchange and correlation effects 
of the FQHE are very well represented by the LDA and that our approach provides 
a computational scheme to model large inhomogeneous FQHE systems. We note that 
there exist previous formal DFTs for strongly correlated systems, in particular for 
high-temperature superconductors [34], and DFT calculations of high-T c materials [35] 
and transition- metal oxides [36]. However, ours are, to the best of our knowledge, 
the first practical LDA-DFT calculations of a strongly correlated system in strong 
magnetic fields, and demonstrate the usefulness of the LDA-DFT in studying large 
inhomogeneous FQHE systems. Recently, we have generalized our DFT approach to 
include spin degrees of freedom in an approximation in which the spin-quantization axis 
is parallel to the external magnetic field. While such an approximation cannot capture 
skyrmion-like excitations, in which the spin quantization axis is tumbling in space, it 
does demonstrate the existence of spin structures near edges of inhomogeneous systems 
consistent with numerical diagonalizations[25] and Hartree-Fock calculations[15]. We 
are presently working on extending the DFT approach to allow for the possibility of a 
tumbling spin quantization axis. 

Ensemble density functional theory approach for spin-polarized systems 

In typical DFT calculations of systems of N e i electrons, the standard Kohn-Sham 
(KS) scheme[37] is implemented, in which the particle density n(r) is expressed in 
terms of a Slater determinant of iV > N c \ KS orbitals, ip a {v). These obey an effective 
single-particle Schrodinger equation H c gip a = e a ip a , which is solved self-consistently 
by occupying the N e \ KS orbitals with the lowest eigenvalues e a , and iterating. This 
scheme works well in practice for pure-state f-representable systems, for which the true 
electron density can be represented by a single Slater determinant of single-particle 
wavef unctions. However, when the KS orbitals are degenerate at the Fermi energy 
(which we identify with the largest e a of the occupied orbitals) there is an ambiguity 
in how to occupy these degenerate orbitals. There exists an extension of DFT which is 
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formally able to deal with this situation. This extension is called ensemble DFT[27, 28], 
and in it, the density of the system is represented by an ensemble of Slater determinants 
of KS orbitals. However, while it can be shown using ensemble DFT that such a 
representation of the density is rigorous, it cannot be shown how the degenerate KS 
orbitals at the Fermi energy should be occupied, i.e., there has not been available a 
practical computational scheme for ensemble density functional theory. 

We begin by first demonstrating that fractional quantum Hall systems are not in 
general pure-state v -represent able, which means their densities cannot be obtained by 
a single Slater determinant of Kohn-Sham orbitals. Hence, an ensemble representation 
has to be used. Consider a uniform fractional quantum Hall system at filling factor 
v = 1/3. This density can be obtained by forming a density matrix of the form 

D=\j:\%)(%\. (2) 

6 i=l 

Here, are the three possible degenerate Slater determinants obtained by filling 
every third single-particle orbital in the lowest Landau level. Using momentum eigen- 
functions, we can write in occupation-number representation = {100100100. . .}, 
|^ 2 > = {010010010...}, and |^ 3 > = {001001001...}. The corresponding ensemble 
density is then 

M r)=TY{i>ft(r)} = ig B ,(r) = 57 i^, (3) 

that is, the filling factor is fixed at v = 1/3. Because we can construct the correct 
ground state density from a density matrix of the form Eq. @ with q > 2, it follows 
from a theorem by Levy [38] and Lieb[39] that this system then has a density which 
cannot be derived from a single ground state Kohn-Sham determinant, i.e.,, no(r) is 
not pure-state w-representable. (This is true whenever we can write the ground state 
density of a system in the form of Eq. (fj) with more than two terms.) However, no(r) 
is still associated with an external potential, and is ensemble representable. 

Although ensemble DFT has been developed formally, there are in practice few 
examples of applications and calculations using ensemble DFT for ground state calcu- 
lations. A significant aspect of our work is that we have developed an ensemble scheme 
which is practical and useful for the study of the FQHE. In ensemble DFT, any physical 
density n(r) can be represented by n(r) = J2 mn fmn\4>mn(r)\ 2 , where f mn are occupation 
numbers satisfying < f mn < 1, and the orbitals ip mn satisfy the equation 



2m* 



e 12 
p + - A(r) 

c 



+ Vext(r)+V H (r) + Kc(r,B)U 

mn(^") CmnV , mn(r), (4) 



where V x A(r) = B(r). In equation @, Vn(r) is the Hartree interaction of the 2D 
electrons, and, as usual, I4 c (r, B) is the exchange-correlation potential, defined as a 
functional derivative of the exchange-correlation energy E xc [n(r), B] of the system with 
respect to density: V xc (r, B) = ^gjfff ■ will hereafter not explicitly indicate the 
parametric dependence of V xc and E xc on B.) For the case of the FQHE, we know that 
the exchange-correlation potential will be crucial, as it contains all the effects of the 
electron correlations which cause the FQHE in the first place, and a major part of the 
DFT application is to come up with an accurate model of E xc and so of V xc . Leaving 
this question aside for a moment, and assuming that we have succeeded in doing so, the 
practical question is then how to determine the KS orbitals and their occupancies in the 
presence of degeneracies. We have devised an empirical scheme, which produces a set 
of occupancies for the KS orbitals which satisfy some minimum requirements, namely 
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(a) the scheme converges to physical densities (to the best of our knowledge) for FQHE 
systems, (b) it reproduces finite temperature DFT distributions at finite temperatures, 
and (c) it reproduces the standard Kohn-Sham scheme for systems whose densities can 
be represented by a single Slater determinant. 

In our scheme, we start with input occupancies and single-particle orbitals and 
iterate the system N cq times using the KS scheme. The number N eq is chosen large 
enough (about 20-40 in practical calculations) that the density is close to the final 
density after N eq iterations. If the density of the system could be represented by 
a single Slater determinant of the KS orbitals, we would now essentially be done. 
However, in this system there are now in general many degenerate or near-degenerate 
orbitals at the Fermi energy. After each iteration, the Kohn-Sham scheme chooses 
to occupy the N e \ orbitals with the lowest eigenvalues, corresponding to making a 
distinct Slater determinant of these orbitals. But there will be small fluctuations in the 
density between each iteration, which cause a different subset of these (near) degenerate 
orbitals to be occupied after each iteration. This corresponds to constructing different 
Slater determinants after each iteration, and the occupation numbers f mn of these 
orbital are zero or unity more or less at random after each iteration. This means 
that the computations will never converge. However, the average occupancies, i.e., the 
occupancies averaged over many iterations, become well defined and approach a definite 
value, e.g., 1/3 for orbitals localized in a region where the local filling factor is close 
to v = 1/3. Therefore, we use these average occupancies to construct an ensemble by 
accumulating running average occupancies (f mn ) after the initial N eq iterations 



where f mn ,i is the occupation number (0 or 1) of orbital ip mn after the ith iteration, 
and use these to calculate densities. Thus, our algorithm essentially picks a different 
(near) degenerate Slater determinant after each iteration, and these determinants are 
all weighted equally in the ensemble. It is clear that this scheme reduces to the KS 
scheme for which the density can be represented by a single Slater determinant of KS 
orbitals (for which the KS scheme picks only the one Slater determinant which gives 
the ground state density) for N eq large enough. Moreover, we have numerically verified 
that a finite-temperature version of our scheme converges to a thermal ensemble at 
finite temperatures down to temperatures of the order of l^huo c /k B . We have also 
performed some Monte Carlo simulations about the ensemble obtained by our scheme. 
In these simulations, we used a Metropolis algorithm to randomly change the occupation 
numbers about our converged solution, keeping the chemical potential fixed. The free 
energy of the new set of occupation numbers was calculated self-consistently. If the free 
energy decreased, this set was kept, and if the free energy increased, the set was kept 
if a random number was smaller than exp [—AF/k B T*], where AF is the change in 
free energy, and T* a fictitious temperature. The results were that to within numerical 
accuracy our ensemble DFT scheme gives the lowest free energy. As a condition for 
convergence, we typically demanded that the difference between the input and output 
ensemble densities, n in (r) and n out (r), of one iteration should satisfy 



Practical density functional theory calculations hinge on the availability of good 
approximations for the exchange-correlation potential V xc , which enters in the effective 
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Schrodinger equation for the KS orbitals. The simplest, and probably the most com- 
monly used, approximation is the local density approximation (LDA). In this approx- 
imation, the exchange-correlation energy is assumed to be a local function of density, 
so that the total exchange-correlation energy consists of contributions from the local 
density of the system. Thus, in this approximation one writes E xc /N = J dre yLC {v)n{r) , 
where e xc (is) is the exchange-correlation energy per particle in a homogeneous system 
of constant density n = v/(2it£ B ) and filling factor v. In other words, in the LDA 
one assumes that the system is locally homogeneous, i.e., the system can locally be 
approximated to have the energy per particle of an infinite, homogeneous system of the 
local density. This approximation obviously makes sense if the density of the system 
varies on a very long length scale, while it could be questionable for systems in which 
the density varies on some microscopic length scale. However, experience has shown 
that the LDA often works surprisingly well, even for systems in which the electron 
density is strongly inhomogeneous[26]. In the FQHE, the length scale of exchange- 
correlation interactions and density fluctuations is given by the magnetic length £ B due 
to the Gaussian fall-off of any single-particle basis in which the interacting Hamilto- 
nian is expanded. The densities are relatively smooth on this length scale, which gives 
us additional hope that the LDA will work well for the FQHE, too. In addition, the 
cusps in the exchange-correlation energy will suppress density fluctuations, so in this 
sense one can actually expect the basic physics of the FQHE to make the LDA a good 
approximation. 

We construct our exchange-correlation energy by writing 

e X cH=erW+ fx C cW. (7) 

Here, e^ M (u) is a smooth interpolation formula (due to Levesque, Weiss and Mac- 
Donald[40]) between ground state energies at some rational fillings. The second term, 
e^cfV), is all-important for the study of the FQHE. This term contains the cusps in the 
ground state energy which cause the FQHE. Here we have used a simple model which 
captures the essential physics. We model e xc (u) by constructing it to be zero at values 
of v — p/q which display the FQHE. Near v — p/q, e^ c {v) is linear and has at v = p/q 
a discontinuity in the slope related to the chemical potential gap A/z = <?(|A p | + \A h \). 
Here A Pi h are the quasiparticle (hole) creation energies which can be obtained from 
the literature [41, 42] at fractions v — p/q. Farther away from v — p/q, e xc {v) decays 
to zero. Finally, in the LDA V xc {r) is obtained from e xc {v) as V xc (r) = ^iai^Ml ^ 

at constant B. In our calculations, we restrict ourselves to include only the cusps at 
v = 1/3, 2/5, 3/5 and v = 2/3 (and the analogous ones at 4/3, 7/5, 8,5 and 5/3), which 
are the strongest fractions. These are some of the fractions of the form v = ( 2p ±i) 
generated by the so-called Vi-model, in which only the pseudo-potential V\ is included. 

A technical difficulty arises in the LDA: the discontinuities in V xc (r) in the LDA 
give rise to a numerical instability. The reason is that an arbitrarily small fluctuation 
in charge density close to an FQHE fraction gives rise to a finite change in energy. 
Imagine that the local filling factor u(r) in some neighborhood of a point r is very 
close to, but less than, say, 1/3 after one iteration. In this neighborhood, the local 
exchange-correlation potential will then form a potential well with sharp barriers at 
the points around r where u(r) = 1/3. During the next iteration, charge will then be 
poured into this well. As a result, the local filling factor will after this iteration exceed 
1/3, and in this neighborhood V xc now forms a potential barrier of finite height. So 
in the next iteration, charge is removed from this neighborhood, and so on. We can 
see that this leads to serious convergence problems. To overcome this, we made the 
compressibility of the system finite, but very small, corresponding to a finite, but very 
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Figure 2. Exchange-correlation potentials (solid line) and V xc i (dashed line) for a 
spin-polarized system as function of filling factor in units of e 2 /(eo^s) for < v < 1.5. The increase 
in V xc at an FQHE filling factoroccurs over a range of filling factor of 0.002. 

large, curvature instead of a point-like cusp in e xc at the FQHE fractions. In other 
words, instead of having a step-like discontinuity A/z in the chemical potential, it rises 
smoothly an amount Afi over an interval 7 in the filling factor. What we found worked 
very well in practice was to have the discontinuity in chemical potential occur over an 
interval of filling factor 7 of magnitude 10 -3 . This corresponds to a sound velocity 
of about 10 6 m/s in the electron gas, which is three orders of magnitude larger than 
the Fermi velocity of a 2D electron gas at densities typical for the FQHE. In general, 
the finite compressibility does not lead to any spurious physical effects so long as the 
energy of density fluctuations on a size of the order of the systems size is larger than any 
other relevant energy in the problem. The only noticeable effect is that incompressible 
plateaus, at which the density would be perfectly constant were the compressibility zero, 
will have density fluctuations on a scale of 7. Figure ^ depicts V xc for a spin-poalized 
system used in our calculations as a function of filling factor. 

Applications to spin-polarized quantum dots 

We have self-consistently solved the KS equations Eqs. @ for a spin-polarized 
quantum dot in a parabolic external potential, V e ^{r) = |m*f2 2 r 2 , by expanding the 

KS orbitals ip mn (r) = e im< ^Lp mn (r) in the eigenstates of H = (p + - c A(r)) . We 
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use the cylindrical gauge, A(r) = |-Br0, and include the four lowest Landau levels 
(n — 0, ... ,3). We chose the static dielectric constant eo = 12.4, appropriate for GaAs, 
and a confining potential of strength[8] HQ — 1.6 meV. 

The use of our LDA-DFT scheme is illustrated by a study of the edge reconstruction 
of the quantum dot as a function of magnetic field strength. As is known from Hartree- 
Fock and exact diagonalizations[21, 22, 25, 18, 19], for strong confinement the quantum 
dot forms a maximum density droplet in which the density is uniform at v = 1 in the 
interior, and falls off rapidly to zero at r « ^2N e \iB = r . As the magnetic field strength 
increases, a "lump" of density breaks off, leaving a "hole" or deficit at about r = r . 
This effect is due to the short-ranged attractive exchange interaction: it is energetically 
favorable to have a lump of density break off so that the system can take advantage of 
the exchange energy in the lump. As B is further increased, the correlations will cause 
incompressible strips with densities v = p/q to appear [23, 24, 43, 31] on the edges, 
and incompressible droplets to form in the bulk at densities u = p/q. Figure [3] depicts 
various stages of edge reconstruction obtained by us as the magnetic field strength 
is increased. The value of B for which the exchange lump appears compares very 
well with the value found by De Chamon and Wen [22] in Hartree-Fock and numerical 
diagonalizations. At higher fields still, incompressible strips appear at the edges, and 
incompressible droplets are formed in the bulk. 

Figure |] depicts the eigenvalues of the KS orbitals for N e \ = 40, and B = 4.10 T. 
The dashed line indicates the chemical potential of the system. This figure then shows 
that all KS orbitals in the bulk are in fact degenerate. It may at first seem paradoxical 
that the eigenvalues are degenerate on an incompressible strip, since, according to 
the picture by Chklovskii, Shklovskii, and Glazman[24], on such a strip the density is 
constant, while the total potential varies (since the electrons cannot screen the external 
potential). If the total potential varies, then ought not the the eigenvalues of the KS 
orbitals localized on that strip vary, too, since these then in general are subjected to 
different potential energies? The problem with this argument as applied to DFT is that 
it ignores the effect of the exchange- correlation potential. As the external and Hartree 
potentials vary across the strip, the exchange- correlation potential varies across its 
discontinuity so as to completely screen out the external and Hartree potentials. The 
discontinuity in V xc does not mean that this potential is fixed at the lower limit of its 
discontinuity while the density is fixed at an incompressible strip. What it does mean, 
is that V xc is free to achieve any value across its discontinuity so as to completely screen 
out the external and Hartree potentials. In this way, it is perhaps better to think of 
incompressibility as the limit of a finite compressibility approaching zero. A strip can 
then remain incompressible with constant density so long as V xc can screen the external 
and Hartree potentials, so the width of the incompressible strip is given by the distance 
over which the external plus Hartree potentials varies an amount given by the energy 
gap associated with the density at that strip. Also, all bulk KS states are degenerate at 
the chemical potential. When a single particle is added, the chemical potential simply 
increases a small amount, and all KS orbitals are again degenerate at the chemical 
potential. 

We have also tested the accuracy of our ensemble DFT-LDA approach by compar- 
ing our results for a six-electron system in a confing parabolic potential with the nu- 
merical diagonalization results by Yang, MacDonald and Johnson[25]. Figure |5| depicts 
angular momentum vs. magnetic field strength for this system. For better compari- 
son with the numerical diagonalizations, we used here only basis states int the lowest 
Landau level (n = 0). There are clear plateau structures in the angular momenta, and 
readily identifiable transitions, such as the initial instability of the maximum density 
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Figure 3. Edge reconstruction of a quantum dot as the magnetic field strength is increased. Plotted 
here is the local filling factor v(r) for a parabolic quantum dot with Ml = 1.6 meV and 40 electrons. 
For magnetic field strengths B < 2.5 T the dot forms a maximum density droplet, and for B w 2.8 T, 
an exchange hole is formed. For stronger magnetic fields, incompressible regions form, separated by 
compressible strips. 



droplet at about 2.8 T, and a formation of a v = 1/3 droplet at about 5.3 T. The agree- 
ment is in general very good, in particular considering that (a) the DFT calculations are 
not constructed to give quantized angular momentum, (b) there are no adjustable pa- 
rameters in our approach, and (c) the numerical diagonalizations used the full Coulomb 
interaction, while our DFT approach only included the cusps at v — 1/3,2/5,3/5, 
and v = 2/3. The Levesque-Weiss-MacDonald interpolation formula used here tends 
to overestimate the magnitude of the exchange-correlation potential around v = 1/2, 
which increases the values at which transitions occur. For example, the v = 1/3 droplet 
formation occurs at 5.29 T in the numerical diagonalization, but at about 5.5T in the 
DFT-LDA calculation. The overestimation of the exchange-correlation potential also 
leads to decreased angular momenta at the different plateaus. For example, while the 
v — 1/3 droplet has angular momentum 45 in the numerical diagonalization, exami- 
nation of the density profiles shows that the formation of a v — 1/3 droplet occurs at 
an angular momentum of about 40 in our calculations. We are presently working on 
improving the exchange- correlation energy for better agreement with numerical diag- 
onalization. Initial calculations give highly accurate results for the maximum density 
droplet instability, and the formation of a v — 1/3 droplet. 
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Figure 4. Eigenvalues of the lowest-Landau level Kohn-Sham orbitals for iV i = 40 and B = 4.10 T 
as a function of angular momentum quantum number. The chemical potential is indicated by the 
dashed line. 



Including spin degrees of freedom 

As mentioned above, due to the small effective Lande ^-factor in GaAs, the spin 
degree of freedom play an important role in quantum Hall systems, in spite of the large 
applied magnetic fields. We have begun to generalize our ensemble density functional 
approach to include the spin degree of freedom. In general [44], spin density functional 
theory has to be based on the single-particle density matrix Pcr CT '(r) = (ipl(r)ip a ' (r)) , 
where ^(r) is the annihilation operator for an electron of spin a at position r. However, 
in the presence of a uniform external magnetic field B = Bz, the z-component of total 
electron spin, S z commutes with the Hamiltonian, and it is a reasonable approximation 
to take Pcrcr'(r) to be diagonal in the spin indices, p CTCr '(r) = n a (r)5 aa i, with n CT (r) the 
up- and down-spin densities. (Note, however, that this restriction has to be lifted in 
order to study skyrmion-like spin textures, which complicates the formalism a great 
deal.) We now obtain two sets of KS equations, one for each spin direction: 



1 



[ 2m* 



p + -A(r) 

c 



+ V; xt (r) + V K (r) + V xc , a (r, B) + ag*fi B Vw(r) 



(r). 



(8) 
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Figure 5. Angular momentum vs. magnetic field strength for a six-electron droplet in a parabolic 
potential with fiQ. = 2.0 meV. The solid line is from numerical diagonalizations (Rcf. 25), and the 
diamonds arc from our DFT-LDA calculations. 



Here, a — ±1, /jl is the Bohr mangeton, and in the local spin density approximation 
(LSDA) the exchange-correlation potentials are 



d 

K c , ff (r, B) = — [ne xc (n T , n b B)] 



(9) 

n rT =n (T (r) 



It is more convenient to work with filling factor u(r) = 2ir£ 2 B n(r) and spin polarization 
£ = (n-f — 7ijJ/(rif + ni), in terms of which we have 

d d 

VW = ^e xc )-(l + 0^e xc , (10) 

where now the exchange-correlation energy per particle in a homogeneous system with 
a filling factor v and polarization £ has to be approximated. Except for a few data 
points obtained by small system numerical diagonalizations [45], this is largely unkown. 
In order to obtain a useful approximation, we start by considering only the exchange 
energy E x \y^, Since the exchange interaction only couples electrons with parallel 
spins, we have 

E xWh v i\ = o#xh, ^t] + a E xWi, nl- (n) 
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In two dimensions, the exchange energy scales as n 3 / 2 , so using Eq. (|TT|), we can 
correctly interpolate between a fully polarized systems (£ = 1) and a completely un- 
polarized one (£ = 0) by writing 

e x [u,£] = e x {y,i= 1) + [e x (u,t = 0) - e x (y,i= 1)] /(£), (12) 

where 

= (l + 0^ + (l-Q^-2^ 
v; 2-2^/2 

We then use the same interpolation for the correlation energy (excluding the cusps for 
a moment), as was first done by von Barth and Hedin[44] for the three-dimensional 
electron gas in zero magnetic field, and write 

4(">£) = <g™(y,t = i) + M",£ = o) - W,Z = i)] /(0 

^6^>,e = i)+5e x »/(e). (14) 

The function 5e xc (z/) can the be obtained by calculating the energy difference between 
polarized and un-polarized systems using data obtained from small system numerical 
diagonalizations [45]. Accurate value for the exchange-correlation energy for £ = and 
v — 1 is not available, but it is a useful approximation to assume that the different spin 
directions are completely uncorrelated at £ = and v — 1, which gives e xc (z/ = 1,£ = 
0) = e xc (z/ = 1/2, £ = 0). This approximation also gives a good value for the exchange 
enhancement, which is the energy required to flip a spin in a polarized system at v = 1. 

So far, we have constructed a function e' KC (v, £) which gives a smooth interpolation 
for the exchange-correlation energy for any value of z/ and £. What is left is to add the 
cusps to this function. We already have a good approximation for this at £ = 1. We now 
need to extend this approximation to arbitrary values of £. Very little is known about 
the cusps, i.e., the energy gaps, for arbitrary polarizations. It is known that there is a 
gap for un-polarized systems at fillings v = 2/5, v = 3/5, and v = 2/3. The gap, and 
thus the cusps, occur at very special 'magic' configurations at which the system can 
take advantage of a particularly low correlation energy. Therefore, it seems plausible 
that for a given value of v, say v = 2/5, there cannot be an energy gap for any value of £ 
between and 1. In order to incorporate this assumption into a usable approximation, 
we interpolate our cusp energy constructed for polarized systems, e xc (^), to arbitrary 
polarizations by multiplying it by a function g(£) which is unity at £ = and £ = 1, 
and vanishes away from these values of polarization. All together, then, we have 

e»(^0 = + + (is) 

Figure depicts V xca (here for £ = 1) used in our calculations. 

We have applied this spin ensemble DFT to study the phase diagram of a maximum 
density droplet. For large values of the Lande-factor g, the maximum density droplet is 
fully polarized, and as the magnetic field is increased, there is an instability to forming 
a spin-polarized exchange-hole. But for small values of g, the instability is towards 
forming a spin structure at the edge. The value of g separating the spin-polarized 
and spin-structured instabilities, g = gfisB/ (e 2 /£q£b) ~ 0.05, is in good agreement 
with the value found for g w 0.03 found from numerical diagonalizations of parabolic 
dots by Yang, MacDonald and Johnson [25], and consistent with the value obtained in 
calculations by Kivelson et al. [15], who used a Hartree-Fock approximation in which 
the spin axis was allowed to tumble. They obtained a value g = 0.17 for an infinite 
Hall bar. 
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Conclusion and summary 

In conclusion, we have showed that ensemble density functional theory can be 
applied to the FQHE. This opens the door to doing realistic calculations for large 
systems. We believe that our results are also significant in that they are the first 
LDA-DFT calculations of a strongly correlated system in a strong magnetic field, and 
they are (to the best of our knowledge) the first practical ensemble DFT calculations. 
We find excellent agreement between our ensemble DFT calculations and numerical 
diagonalizations and Hartree-Fock calculations. Preliminary calculations including spin 
degree of freedoms are consistent with numerical diagonalizations and Hartree-Fock 
calculations. We are presently working on improving the exchange-correlation energy 
as a function of density and polarization. Finally, the spin DFT presented here cannot 
be used to study skyrmions-like structures, in which the spin density vector is tumbling 
in space. Work is presently underway, together with J. Kinaret (Chalmers University 
of Technology) to develop such a theory. 
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O.H. would like to thank Chalmers Institute of Technology, where part of the numerical 
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